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Abstract. We present a A''-diinensional generalization of the two-dimensional block-circulant perfect 
array construction by [Blake, 2013| . As in |Blake, 2013] , the families of A'^-dimensional arrays possess 
pairwise good zero correlation zone (TjCL) cross-correlation. Both constructions use a perfect autocorre- 
lation sequence with the array orthogonality property (AOP). 

This paper presents a generalization of the 2-diinensional block-circulant array construction by [Blake, 2013| . 
. . , The generalized construction works in any number of dimensions, but is limited to the same size in each 

CO ■ dimension as the original two-dimensional construction. Throughout the paper we follow the notation of 

'^ \ [Blake, 2013| . 

(N . 

, We begin by stating the 2-dimensional construction for families of arrays with perfect autocorrelation 

O ■ and good cross-correlation as given in [Blake, 2013| . 

Construction I Let a — [ao,ai, • • ■ ,a„_i] be a perfect sequence with the AOP for the divisor, d, and 
'^^ ' and c = [c(0), c(l), • • ■ , c{d — 1)] is a block of d perfect sequences - each of length m, where m = mod d. 

We construct a family of arrays, Sfc, such that 

Sfe = [Si,j]k ~ 0,j c{j mod d).ii,Lj7dJ+A;(j mod d)+i 

Q ■ for < i < n, < j < m, a has the AOP for the divisor d, < k < m, and w = m/d. 

This construction produces perfect arrays up to size r^ x r^ over r roots of unity. Each pair of distinct 
^ ] arrays has d^ non-zero cross-correlation values, as d <C r we say the array has good ZCZ cross-correlation. 

OO ' The generalized A^-dimensional construction is given as follows. 

0\ 

Construction II Let a = [ao,ai, ■ ■ ■ ,a„_i] be a perfect sequence with the AOP for the divisor, d, and 
JnI i c = [c(0), c(l), • • • , c{d — 1)] is a block of d perfect sequences - each of length m, where m = mod d. We 

construct a family of k A^-dimensional perfect arrays, Sfc, such that 



o 



X 



N-2 

Sfc = [5*40,41,- ,iN_2j]fc = "-j 11 c(j mod d).u,ij/d]+k{j mod d)+j„ 

11=0 

for < io,ii-, ■ ■ ■ ,i7V-2 < n, < j < m, < k < m, and w — m/d. 

Each distinct pair of arrays from the family has d^ non-zero cross-correlation values, regardless of the 
size of N. Thus, the ratio of zero to non-zero cross-correlation values is larger for higher dimensional arrays. 

We show that Sfc from Construction II has perfect autocorrelation (fci = fc2) and S^^, S^^ has good 
cross-correlation (fci ^ ^2). We now compute the cross-correlation of S^-^ with Sfc^, (0 < ki,k2 < m) for 
shift so,si, • • • ,SN-V 



io + so,ii+si,--- j'+sjv-l 



9Sfc,,Sfc, (so,.Si,--- ,SAr_i) 
Ti— 1 ni—1 ?7i— 1 rn — 1 

ri-lm-lm-1 m-1 / /N-2 \ 

= XI 51 H ■ ■ ■ H I aj I n '^(■^' ™°^ C?)^Lj7dJ+fci(j mod d)+i„ j aj+,,„_i X 

( 11 C{j + SN~1 mod rf)Jl,L(j + ,,„_i)/rfJ+fc2(i+SN-l modrf)+j„+.s„ 1 1 

n-1 / m-1 m-1 )ri-l / /JV-2 \ 

= H «J aj + s«_i H II ■ • • H n ^(^' "^°^ rf)«.Lj7<iJ+fel(j mod d)+^^ X 

i=0 y io=Oii=0 ijv_2=0 \ \i;=0 / 

n c(j + SAr-l mod rf);L(j + s„_i)/dJ+fe2(j + s«-l modd)+.„+«„ 



(7) 



Consider the case when sn-i ~ mod d. Let sjv-i = I d and perform the change of coordinates j = qd + r, 
(r < d). Then the inner summation in (7) becomes 



m— Im— 1 m — 1 / /N~2 



/JV-2 



E E • ■ • E n ^('^)^-+-^ n ^(0L+.2.+.„+.s„ 

io=Oii=0 i„_2=0\\t)=0 / \?)=0 // 

(m-1 \ /m-1 \ 

E c('')fcir+»o c{r)*i^+k2r+to+so ) f E c(^)fcir+ii c{r)L+h2r+t,+.n ) ^ 
20=0 / \n=0 / 

/m— 1 \ / m— 1 

I 2_^ c(r)/j;jr+i2 c(r);^^j,2^_|.j^_,_j^ 1 X • • • X I 2_^ c{r)kj^r+iN-2 c(?');„,+fe2r+'ijv-2+s 



\i2=0 



\iN^2=0 



which is independent of q. Thus, (7) becomes 

d-l / Ai/d-l \ /m-1 \ /m-1 ^ 

r=0 V V 9=0 / \io=0 



^ii=0 



/m— 1 



m—1 



/ , c(r)fejr+i2 c(»');„,+fe2»-+ 



'i2+S2 



7 ^ C(rjfcj^r+iiV-2 C(?');-uj + fc2'- + iN-2+S 



Vi2=0 



i.iiv-2=0 



(8) 



As a has the AOP for the divisor d, the summation X]g=o "'3<i+'' "qd+r+sw-i = ^ f^"" SAr_i ^ mod n/d, 
thus 0Sfc ,Sfc (sojSii"'- ,SAf-i) = for sat-i 7^ mod njd. Otherwise, for s^v-i = mod n/d, (8) is zero 
when one of the remaining N ~1 summations are zero. The summation over i^ is non-zero for (ki — A:i)r + 
sjv-i'Ti/d^ + s„ = mod m. For each of the sjv_i = 0, n/d, 2n/d^ ■ ■ ■ , (d — l)n/d, there are d solutions for 
Su, one for each r = 0, 1, • • • , d— 1. Thus, there are d^ non-zero cross-correlation values for Sat-i = mod d. 



Now consider the case fci = fc2 (autocorrelation), (note that sn-i = mod n/d), then (8) becomes 

^d-ln/d-l \ /,n-l \ /,ri-l \ 

^r=0 g=0 y \io=0 / \»i=0 / 

C)ri-1 \ / m-1 \ 

J2 ^W»^ (=ir)L+^,+s, P • • • ^ II c(r),„_, c(r);„+,„_^+,^_^ (8) 

As a has the AOP for the divisor rf, the left double summation is zero. Thus the autocorrelation is zero for 
Sjv-i = mod d. 

Now, consider the case when sjy-i ^ mod d. Let sn-i = I d + s, (s < d), and perform the change of 
coordinates j = qd + r, (r < d). Then (7) becomes 

d-m/rf-l / m-im-i m-1 / /N-2 \ 

r=0 g=0 y io=OJi=0 ijv-2=0 \ \ u=0 / 

I 11 c(r + S mod rf)q„,+;„,+^L(r+is)/dJ+fe2(r+s mod d)+j„+s„ I I I 

d-ln/d-1 I /m-1 \ 

= /_^ Z_^ I '^ijd+r aqd+r+sjv-i I 2^ c(r)q.„+fcjr+io '^V + ^ l^^od ")gu;+Zu)+«) [(■'-+s)/dj +'i:2('-+s mod d)+io+so I ^ 
l-=0 q=0 \ \lo=0 J 

I 2_^ c(r),u,+fc^r+ii c{r + S mod rf)gu,_|_/„+^L(»-+s)/dJ+*;2(i-+s mod d)+ii+si I X ■ ■ • X 
\ii=0 / 

'""^ * ^ ^ 

/ , '^(^)9"'+'i:i''+«jv-2 c(?' + s mod d)^^_,_;^^_,^,^(^^^j/^j^j,^(,._i_^ jj^qj ^^^j^_^_^j^_2 I I (9) 

JV-2=0 J J 

For each of the inner summations above, shifting both sequences qw places does not change the cross- 
correlation. So (9) becomes 



d-l I /n/rf-1 

//-(iJ-f -u; L{r+s)/rfJ +^2 (^'+s mod d)+'io + SG 



r=0 \ \ 9=0 ) \Jo=0 

/""^ * \ 

I Z^ c(r)fejr+ii c{r + S mod rf);*u,+„,|^(r+s)/dj+fe2(r+s mod d)+ii+si I X ' ' ' X 

\n=0 / 

ra—l 

/ . c(r)A;jr+i„_2 c(r + S mod rf)L+„L(r+s)/dJ+fc2(r+s mod d)+iiv_2+sjv-2 
^ijv-2=0 

As the sequence a has the AOP, the summation 

n/d-l 
<j=0 

for all r and sat-i such that sat-i 7^ mod rf. Thus the cross-correlation of 8/;^ with S^-j is zero for 
h ^ mod d and all z;. So Construction II produces a family of m perfect fc-dimensional arrays with pair- 
wise good cross-correlation. 



It has been shown that perfect iV-dimensional arrays can be constructed up to size r^ x r^ x • • • over 
r roots of unity. An exponential number of such arrays exist for each size. Furthermore, famihes of arrays 
possess good cross-correlation properties. Importantly, despite the size of iV, the number of non-zero cross- 
correlation values is always (P , so the ratio of the number of non-zero cross-correlation values to the number 
of entries in the arrays can be made as small as desired. 
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Appendix I — Implementation of the Construction 

In this appendix we give a IMathematical (version 8.0) implementation of the construction. (Arrays are given 
in index notation, that is, the mapping: Q'^'^y-^sr^/r _^ ^^^ -j 

We begin with the code for periodic cross-correlation, XCV and autocorrelation, ACT: 

In[l]:= XCV[a_, b_, r_Integer] := Block [-[A, B}-, 
A = Developer 'ToPackedArray [Exp [(2 . Pi I a)/r]]; 
B = Developer'ToPackedArray[Exp[(-2. Pi I b)/r]]; 
Chop[ListCorrelate[A, B, 1], l*~-5]] 

In[2] := ACV[m_, r_Integer] := XCV[m, m, r] 

The A^-dimensional array construction, Construct ionND, takes as input the multiplying sequence, a, the 
block of d perfect sequences, c, the parameter, fc, the number of dimensions, dims, and the number of roots 
of unity 

In[3]:= aref = #1[[#2 + 1]] ft; 

In [4]:= ConstructionND [a_, c_, k_, dims_] : = 
With[{d = Length [c], 1 = Length[First [c]] >, 
Table [ 
Array [ 

indexND[a, c, k, d, 1, j, ##] &, Table[l, {dims - I}-], 0], 
{j, 0, Length [a] - 1}]] 

In [5] : = indexND [a_ , c_ , k_ , d_ , 1_ , index ] : = 

With[{w = Length [First [c]]/d, j = First {index}}, 
aref [a, j] + Product [aref [aref [c , Mod[j, d] ] , 

Mod[w Floor [j/d] + k Mod[j, d] + i, 1]], {i, Rest {index}}]] 

We construct a perfect 4x4x4x4 binary array and compute its periodic autocorrelation: 

In[6] := ConstructionND [FrankSequence [2] , {#, Decimate [#, 3]} & FrankSequence [2] , 0, 4]; 

Mod[y., 2] 

Count [ACV [•/. , 2], n_?NumericQ /; n != 0., Infinity] 



Out [7] = 
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Out [8]= 1 

The following generalization of indexND allows for computations over the quaternions: 

In [9]:= NonCommutativeDot [e_] := e 

In [10]:= indexND[a_, c_, k_, d_, 1_, index ] /; ! FreeQ[a, Quaternion] : = 

With[{w = Length [First [c]]/d, j = First {index}}, 
aref[a, j] ** Apply [NonCommutativeDot , 

Table [aref [aref [c , Mod[j, d] ] , Mod[w Floor [j/d] + k Mod[j, d] + i, 1]], i±, Rest@-[index}}] ] ] 

The following is a perfect quaternion sequence which posesses the array orthogonality property for the 
divisor 4: 

We use this sequence to construct a perfect 16 x 16 quaternion array. 

In[ll]:= ql6 = {Quaternion [1, 0, 0, 0], Quaternion[0, 0, 0, 1], 
Quaternion[l, 0, 0, 0], Quaternion[0, 0, 0, -1], 
Quaternion [0, -1, 0, 0], Quaternion [0, 0, 0, -1], 
Quaternion[0, 1, 0, 0], Quaternion[0, 0, 0, -1], 
Quaternion [-1, 0, 0, 0], Quaternion[0, 0, 0, 1], 
Quaternion [-1, 0, 0, 0], Quaternion[0, 0, 0, -1], 
Quaternion[0, 1, 0, 0], Quaternion[0, 0, 0, -1], 
Quaternion [0, -1, 0, 0], Quaternion [0, 0, 0, -1]}; 



In[12]:= ConstructionND[ql6, {#, Decimate[#, 3], RotateRight [#, 2], #} & ql6, 0, 2] 
Count [QuaternionCV2D [•/.] , n_?NumericQ /; n != 0., Infinity] 
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Out [13]= 1 



































We now construct a family of 9 perfect 9x9x9x9 arrays over 3 roots of unity. The perfect Frank 
sequence: [0, 0, 0, 0, 1, 2, 0, 2, 1], is used as the base sequence. 



In [14]:= familySD = Table [ 

ConstructionND [FrankSequence [3] , 

■[Decimate [#, 2], Decimate [#, 5], Decimate [#, 7]} 
{n, 9}]; 



FrankSequence [3] , n, 4], 



Each pair of sequences from the family have (P = 3^ non-zero cross-correlation values. Here's the non-zero 
cross-correlation values of an example pair: 

In [15] : = XCV [f amily3D [ [1] ] , f amily3D [ [2] ] , 3] ; 
Cases['/,, n_?NumericQ /; n != 0., Infinity] 



0ut[16]= -[2187., 2187., 2187., -1093.5 - 1894. I, -1093.5 + 1894. I, 2187., -1093.5 + 1894. I, 
-1093.5 - 1894. I, 2187.} 

Here we generate all pairs of arrays from the family and count the number of non-zero cross-correlation 
values for each pair. 



In [17]:= Count [XCV [#1, #2, 3], n_?Numericq /; n != 0., Infinity] 



Tuples [family3D, 2] 



0ut[17]= {1, 9, 9, 9, 9, 9, 9, 9, 9, 9, 1, 9, 9, 9, 9, 9, 9, 9, 9, 9, 1, 9, 9, 9, 9, 9, 9, 9, 9, 
9, 1, 9, 9, 9, 9, 9, 9, 9, 9, 9, 1, 9, 9, 9, 9, 9, 9, 9, 9, 9, 1, 9, 9, 9, 9, 9, 9, 9, 9, 9, 1, 
9, 9, 9, 9, 9, 9, 9, 9, 9, 1, 9, 9, 9, 9, 9, 9, 9, 9, 9, 1} 



(The I's above correspond to the autocorrelation of non distinct pairs.) 



